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ABSTRACT 

We report a Spitzer/IRAC search for infrared excesses around white dwarfs, 
including 14 newly-observed targets and 16 unpublished archived stars. We find a 
substantial infrared excess around two warm white dwarfs - J220934. 84+122336. 5 
and WD 0843+516, the latter apparently being the hottest white dwarf known to 
display a close-in dust disk. Extending previous studies, we find that the fraction 
of white dwarfs with dust disks increases as the star's temperature increases; for 
stars cooler than 10,000 K, even the most heavily polluted ones do not have 
~1000 K dust. There is tentative evidence that the dust disk occurrence is 
correlated with the volatility of the accreted material. In the Appendix, we 
modify a previous analysis to show that Poynting-Robertson drag might play an 
important role in transferring materials from a dust disk into a white dwarf's 
atmosphere. 

Subject headings: circumstellar matter - minor planets - white dwarfs 



1. INTRODUCTION 



White dwarfs are the final evolutionary stage of stars with masses less than about eight 
solar masses. The first infrared excess around a white dwarf was discovered more than 20 



years ago (IZuckerman fc Becklinl 119871 ). Since 2005, primarily thanks to the Spitzer Space 
Telescope ( Werner et al. 2004 ). progress has been dramatic. Including WD 0843+516^ and 
J220934.84+122336.5 (hereafter J2209+1223), which are reported in this paper, there are 23 
confirmed white dwarfs with a warm disk and 5 additional candidates (see our Table [U an 
extension of Table 1 in iFarihi et al.l (120091 ) ). Recently, observations from the WISE mission 
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1 After reading a short description on astro-ph of a survey on dust disks around white dwarfs ( Chu et al 



2010I ). we learned from Y.-H. Chu that their team has independently found the infrared excess for this star. 
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(Debes et al. 


201 J) 


Girven et al. 


2011) 



There are three classes of dust disks detected around a white dwarf: (1) ~100 K dust, 
which is found aroun d > 15% pre-wh ite dwarfs and white dwarfs, sometimes in the center of 



a planetary nebulae (IChu et al.ll201ll ). The source of th e dust is believed to be t he collisions 



(Bonsor & Wvatt 


2010 


: iDone et al. 


I (T*=7,400 K) (I 


? arihi et al. 


2008bl) 



and PG 1225-079 (T,=10,600 K) (IFarihi et al.ll2010b[ ). A model that reproduces these data 
invokes emission from an opaque dust ring with a large, dust-free inner hole and an outer 
boundary within the tidal radius of the white dwarf. However, there are only two dust rings 
in this class and their origin is uncertain. (3) ^1000 K dust, which is found in the remaining 
21 white dwarfs listed in Table [TJ In these systems, all the dust lies within the star's tidal 
radius and the inner boundary is determined by the location where the refractories rapidly 
sublimate. The excess can be fit with a flat opaque dust disk (jJurall2003l ). The excesses from 
1000 K dust are substantial at wavelengths greater than 3 /xm and can be easily studied with 
warm Spitzer, we focus on discussing this class of disk in this papeJE 

Gravitational settling of heavy elements is usually so effective in white dwar fs cooler 
than 25,000 K that their atmospheres are typically pure hydrogen or pure helium (jrCoester 
20091 ). However, about 25% of DAs (white dwarfs with hydrogen dominated atmosphere) and 
33% of DBs (white dwarfs wi th helium dominated atmosphere) have atmospheric pollution 
(jZuckerman et al.ll2003l . |2010| ). Given that the diffusi on timescale o f heavy elements is orders 
of magnitude less than a white dwarf's cooling age (jKoesterll2009i ). it is unlikely that those 
heavy elements are intrinsic. Theoretical calculations show that tidal disruption and subse- 
quent accretio n of small bodies in the planetary system can explain the s o urce of pollution in 
these objects JDebes k SigurdssorJ liool I Jural 120031 : bonsor et aDboilh . iKilic et all (bood ) 
first showed a correlation between the presence of an orbiting dust disk and a high level of 
heavy element pollution in the star's photosphere. Subsequen t observations have confirmed 



this tr en d and suppo rted the small body accretion model (IJura et al.l 120071 ; iFarihi et al. 



200830, 1 20091 . l2010bl ). Thus studying polluted white dwarfs provide s invaluable measu res 



of the bulk composition of extrasolar minor planets (see, for example, iKlein et al.l (120101 ) ). 



However, as can be seen in our Tableland in IFarihi et al.l (120091 ). almost all the disk- 
host stars have effective temperatures of at least 9,500 K. Here, we report the results from 
a warm Spitzer Cycle 7 program, searching for infrared excesses around polluted cool white 



2 In the following sections, unless specifically stated, dust disk and infrared excess all refer to this 1000 K 
dust. 
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dwarfs and a few highly-polluted warmer targets. Additionally, we reduced some unpublished 
archived Spitzer data. 

Target selections and observations are presented in section 2. In section 3, we describe 
data reduction procedures and SED fits to the data. In section 4, the main result is discussed: 
we continue to find dust disks around warm polluted white dwarfs but no dust around highly 
polluted cool stars. In section 5, we explore the possible correlations between the presence 
of a dust disk and other characteristics of the star in addition to its temperature. Possible 
scenarios that can account for the presence of heavy elements without producing a dust disk 
are discussed in section 6. Conclusions are presented in section 7. In the Appendix, we assess 
the importance of Poynting-Robertson (P-R) drag on the entire dust disk as a mechanism 
to transfer material from the disk into the white dwarf's atmosphere. 



OBSERVATIONS 



2.1. Source Selection 



The 14 targets of our Cycle 7 program are listed in Table [2J 10 stars are newly identified 
DZa§ with distinctively large amount of atmospheric calcium det ected in the SPSS spectru m; 



they all have T* < 10,000 K and m(z) < 18.0 mag (0.23 mJy) (IDufour et al 



2006 



2007 L_4 



additional DBZg9 with T* ~ 15,000 K are also included: 3 stars from Eisenstein et al. ( 2006 ) 
plus G 241-6, which was reco gnized recently to be a near-twin of GD 4 0, a highly polluted 
white dwarf with a dust disk (IZuckerman et al.l l2010l : Klein et al.l 120101 ) . 



For the completeness of the sample, we also report IRAC 3.6 /im - 7.9 //m fluxes for 16 
single DA white dwarfs from the unpublished archived Cycle 3 program 30856. These stars 
all have well determined temperatures, surface gravities as well as highly accurate SDSS 
photometry, which is crucial in determining the presence of a disk. Relevant parameters are 
listed in Table [3j 



2.2. Observing Strategies 

Spitzer has been operated in the warm phase mission since the depletion of cryogen 



3 DZs are cool white dwarfs that display trace elements other than carbon; it is thought that DZs have a 
helium dominated atmosphere, but the gas is too cool for the helium lines to be detected. 

4 DBZs are DBs with detected heavy elements in the atmosphere. 
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i n May, 2009. Th e two shortest wavelength channels (3.6 /zm and 4.5 fim) of the IRAC 
( IFazio et al.l 120041 ) continue to function and were used in this study to obtain broadband 
images. The observations were performed by using a 30 sec frame time with 30 medium 
size dithers in the cycling pattern, resulting in a 900 sec total exposure time in each IRAC 
channel. 



The archived data were obtained between 2006 and 2007. The observing mode was 
random 9-point large dithering with a frame time of 30 sec, resulting in a total of 270 sec in 
each IRAC channel. 



3. DATA ANALYSIS 



3.1. Data Reduction and SED Fits 



Following the data reduction procedures described in lFarihi et al.l (j2008al Jbl. 120091 ). each 
exposure was processed with MOPEX (version 18.4.9) to create a single mosaic with a pixel 
size of 0'/6. Though the intrinsic plate scale is 1'.'2 pixel -1 , a 0'.'6 pixel -1 mosaic was used 
because it has a better spatial resolution and the 3.6 [im and 4.5 fim images are oversampled. 
Aperture photometry was performed on the combined mosaic using both: (i) the standard 
IRAF apphot task, and (ii) the Astronomical Point source Extractor (APEX) in MOPEX. 
We tried an aperture size of both 2 (2'.'4) and 3 (3'/6) native pixels with a sky annulus of an 
inner radius of 5 native pixels (6'/0) and an outer radius of 15 native pixels (18'.'0). We found 
that in some cases even in a clean field, the measured flux in these two aperture sizes can 
differ up to 5% and we report the average flux weighted by their SNR. To derive the aperture 
correction factor, we performed aperture photometry in the point response function (PRF) 
images in each band with our set of parameters. The PRF is the convolution of the point 
response function (PSF) of the telescope and the pixel response function of the detector and 
it is provided by the Spitzer Science Center. There are three factors contributing to the 



total u ncertainty: 5% call 



dwarfs (Farihi et al. 



2008a 



Dration uncertainty, which is based on previous studies of white 
), rather th an the 2-3% value optimally derived for other targets 



( IReach et al.l 120051 ; iBohlin et al.l 120 111 ), measurement error and uncertainties caused by the 
choice of aperture radius; these uncertantities were all added in quadrature. The fluxes 
independently obtained by IRAF and APEX agree to within 3%. To be conservative, we 
report the values that have the larger measurement error in Table H] and Table |5j 

When there is a background source, proper motion analysis is performed to accurately 
determine the position of the target, followed by PRF fitting to obtain the true flux of the 
star. A residual mosaic, which subtracted the detected point source from the original mosaic, 



-5 - 



was produced. To check the result of PRF fitting, aperture photometry was re-performed 
at the same location of the target in the residual image. A good fit is achieved if the flux 
level in the aperture is comparable to the sky level, which is a few percent of the flux of the 
source. Examples of PRF fitting and the residual image are shown in Figure [T] and Figure [2J 

For some archived stars, the background is very noisy in the 5.7 /im and 7.9 /im band. 
We report a positive detection when the measured flux is at least 3cr. To measure the upper 
limit, aperture photometry was performed in 20 empty background regions around the source 
with an aperture size of 2 native pixels. Then the standard deviation of these fluxes is taken 
to equal la. 

The Spectral Energy Distribution (SED) of all the targets are plotted in Figures [31 [TOl 
Also included in the SEDs are ugriz fluxes from the SDSS, JHK fluxes from the Two Micron 
All Sky Survey (2MASS) and WISE fluxes wh en the data are ava ilable. Ultraviolet fluxes 
from the Galaxy Evolution Explorer (GALEX) (IMartin et al.ll2005l ) are usually excluded be- 
cause they can be strongly suppressed due to heavy element blanketing (IKoester et al.ll201ll ) 
or interstellar extinction. A blackbody model is then adopted to fit the star's photospheric 
flux, with most weight given to the SDSS photometry because the 2M ASS fluxes have larger 
uncertainties. Though the blackbody temperat ure does not always agree with the reported 
white dwarf temperature of iDufour et al.l (120071 ). this method is sufficient to identify infrared 
excess that is 10% above the photospheric value. When the re is a detected excess, a thin, 
opaque, passive dusty disk is used to fit the SED (I Jural 120031 ) . 



3.2. Notes on Individual Stars 



3.2.1. WD 0216-095 



This star is heavily blended with a nearby galaxy SDSS J021836. 98-091955. 7 with an 
angular separation of 2 / . / 6 at position angle 210°. We ta ke a proper motion of 130 mas/yr 
in RA and 43 mas/yr in Dec from the NOMAD catalog (IZacharias et al.l 120051 ) to derive a 
position of 02:18:36.77-09:19:44.24 at the epoch of the Spitzer observation. Then we used 
the method described above to resolve these two objects. No excess is detected for this star. 



3.2.2. WD 0936+560 



In Figure [3X we see the 2MASS H band upper limit is significantly lower than our 
predicted photospheric flux. Based on the SDSS photometry and the near infrared spectra 



for this star in lKilic et al.l ( 120081 ). we assign a very low weight for this point when fitting the 
SED and no infrared excess is found for this star. 



3.2.3. J2209+1223 

The IRAC image of this star reveals an unresolved background star separated at 6'.'2 at 
position angle 333° (see Figured]). The proper motion of this star is small enough that it 
is neglected. Robust centroid methods were used in each exposure to accurately determine 
the position of these two objects, and then PRF fitting was performed. To demonstrate the 
effectiveness of the PRF fitting, the original image and residue after subtracting the star are 
shown in Figure [TJ 

Figure |9] shows the SED for J2209+1223; we see a strong infrared excess in both IRAC 
ba nds, which is better explained by a flat disk rather than a companion. In the disk models 



of I Jural ( 120031 ). there is a degeneracy between the inclination of the disk and its size. With 
only two data points, there are two sets of parameters as listed in Table |6]that can reproduce 
the SED: a more inclined larger disk or more face-on smaller disk. Another good fit to the 
SED is a blackbody with T = 1000 K and R = 9R at the same distance as the white dwarf, 
190 pc. However, this derived radius is much too large for a brown dwarf. Furthermore, 
there is no feature indicative of M dwarf companions, such as Ha emission or TiO absorption 
lines in the SDSS spectra of this star. 

J2209+1223 has an effective temperature of 17,300 K with an accretion rate of 3 x 
10 10 g s _1 ; this value is quite uncertain as it is derived from scaling the equivalent width 
of magnesium line s from GD 61, which is a white dwarf with similar stellar parameters 



( IFarihi et al.ll201lf ). This star falls into the catalog of highly polluted warm white dwarfs 



with an infrared excess. 



3.2.4. G 241-6 



We take the prop er motion of G 241-6 as 144 mas/yr in RA and 244 mas/yr in Dec 
( iZacharias et al.ll2005l ) to derive its position at 22:23:33.38+68:37:26.77 in the Spitzer obser- 
vation. In Figure |2l we can see in both IRAC bands this star is heavily blended with an 
unknown object at an angular separation of 2'.'3 at position angle 170°. Another galaxy is 
present at J222334. 23+683723.8 at an angular distance of 5'.'5 and position angle 110°. By 
using PRF fitting, we successfully resolve these ob jects and no excess is f ound even though 
this 15,300 K star is as heavily polluted as GD 40 (IZuckerman et al.ll2010l ). which does have 
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a dust disk (IJura et al.ll2007l ) 



3.2.5. WD 0819+363 



In four IRAC bands, this star is blended with the galaxy SDSS J082245.87+361410.2 
with a distance of 3'/5 at position angle 240° and another unknown object separated at 4'.'1 
and pos ition angle 110° . The proper motion of this star is -78 mas/yr in RA and -58 mas/yr 
in Dec (IZacharias et al.ll2005l ). and it is located at 8:22:46.2+36:14:12.04 at the epoch of the 
Spitzer observation. PRF fitting was performed in the 3.6 /im and 4.5 /im band and no excess 
was found. In the 5.7 /im and 7.9 /im band, the background is so noisy that APEX failed to 
perform PRF fitting in the nominal position derived from IRAC band 1 and 2. Instead we 
used aperture photometry at the position of the star and there appears to be a 3cr excess as 
shown in Figure [6j However, this is unlikely to be real considering the complicated field. 



3.2.6. WD 0843+516 

Fig HO] shows the SED for WD 0843+516, apparently the hottest white dwarf so far 
detected to have a dust disk. The model atmosphere was provided by D. Koester (private 
communication). An opaque disk model ca n reproduce the data with the parameters listed in 
Table |6j As previously found with GD 362 ( Jura et al1l2007 ). the model tends to underpredict 
the flux at 7.9 /tm, which might be contaminated by strong silicate emission in t his band. An 



M dw arf companion can be ruled out as there is no excess in the JHK bands (see iHoard et al. 
( 120071 ) ) . Another fit to the SED can be achieved with a blackbody of T = 800 K and R 



O.22R at 80 pc. However, this radius is much too large for a brown dwarf, so this possibility 
is also excluded. So far WD 0843+516 is the only white dwarf with an infrared excess that is 
yet to know whether it is polluted as there is no suitable spectra reported in the literature. 



RESULTS 



4.1. Missing Planetary Debris Around Cool Stars 



Warm dust has been found around 21 white dwarfs with temperatures ranging from 
9,500 K to 24,000 K and at least 20 of them have a highly polluted atmosphere. Between 
1% to 3% single white d warfs with cooling ages less than 0.5 Gyr possess warm circumstellar 
dust ( IFarihi et al.ll2009h . In contrast, the disk fraction around cool white dwarfs (T* < 6000 
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K) only has an uppe r limit of 0.8% (IKilic et al.ll2009l) . Previously, 11 DZs have been targeted 
with Spitzer/IRAC (IFarihi et al.ll2008al ibl. 120091 ). but none shows an infrared excess. In this 
study, we observed 10 additional DZs, which are heavily polluted and most likely disk-host 
stars. Though almost doubling the numbers of DZs, we still fail to find any infrared excess. 
We address the implications and possible explanations in section 5 and 6. 



4.2. Hot White Dwarfs with a Close-in Disk 



As discussed by Ivon Hippel et al.l (120071 ). we need to understand how solid dust grains 
can survive around hot white dwarfs. Dust disks have been confirmed around four stars 
with stellar tempe rature higher t han 20,000 K: WD 0843+516 (this pap er), GALEX 1931 



flDebes et al.ll201lh . PG 1457-086 (IFarihi et al.l 120091 1 and WD 1226+110 (IBrinkworth et al. 



20091 ). Currently, there ar e two differ ent models to describe the inner boundary of the disk: 
a fully opaque disk model (jJurall2003l ) and an opaque disk wi th a thin transition zone, which 
goes from optically thin to optically thick (iRafikovl 12011 al l. However, for the hot white 
dwarfs, only the opaque disk model allows the inner boundary of the disk to be within the 
tidal radius, 



-) 

3tt J 



1/3 



4/3 



where T* and r* is the stellar temperature and radius. T s is the temperature at which the 
dust particles sublim ate; it can be cons trained from the 3.6 ^m excess and it is usually 1200 
K (see, for example, lJura et al.l (120071 )). So for a 20,000 K star, the inner boundary for a 
complete opaque disk is 25 r*. 



In the model described by iRafikovl (j2011al ). the inner boundary of the optically thin 



region is, 



(2) 



_ n (T* 

r inner 2 \ T 

Fo r a 20,000 K star, the in ner radius is 140 r*, whereas the tidal radius is approximately 130 
r* (Ivon Hippel et al.l 120071 ). In this case, all the dust particles lie outside of the tidal radius. 



We have no way to resolve the disk and directly tell which model is correct, but the fully 
opaque disk model is the simplest explanation for the presence of close-in dust. There are two 
additional arguments in favor of the opaque disk model: it can reproduce the substantial 
amount of infrared excess and the dust disk radii agree with those derived from gaseous 
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emission lines (iMelis et al.ll2010l ). Further discussion is in present ed in the A ppendix, where 
we consider P-R drag on a completely opaque disk described by lJural (120031 ) and derive an 
accretion rate t hat agrees better with the data than the model described in Equation (j2j) 
rtRafikovl lioilal l . 



5. CHARACTERISTICS OF DISK-HOST STARS 

We see the planetary debris seems to be missing around the polluted cool stars. Why is 
that? In this section, we explore different characteristics that might be correlated with the 
occurrence of a dust disk. 



5.1. High Heavy Element Accretion Rate? 



It is found that at least 50% of white dwa rfs with heavy element accretion rate over 3 
x 10 8 g s _1 have warm dust (IFarihi et al.ll2009l ). But cool stars do not appear to follow this 
pattern. Figure [IT] shows a comparison of the overall accretion rate with a white dwarf's 
stellar temperature. We see that there are numerous cool white dwarfs (T* < 10,000 K), 
that have an accretion rate higher than 3 x 10 8 g s -1 , but none displays a dust disk. 

To calculate the accretion rate of element A, we assume a steady state so it equals the 
mass of A in the convective zone divided by its settling time. For the stars that have well 
determined major element abundances, the accretion rate is simply the sum of all the heavy 
elements. For the stars that have both magnesium and iron abundance, we assume they 
are 30% of the total mass of the parent body. For the Cycle 7 targets, calcium is the only 
de tected element and w e assume it is 0.6% of the total mass, extrapolated fr om the c o ol DZ s 



in 



Koester et all ( 1201 it ). 



Koester et al.l (120111) . All the stellar parameters are extrapolated from iKoesterl ( 120091 ); 



5.2. Massive Parent Bodies? 



Another possibility is that while more massive parent bodies can form a dusty disk, less 
massive paren t bodie s can only form a more tenuous, short-lived dust disk or a purely gaseous 
disk (see I Jural (120081 )). A comparison of the mass of heavy elements in the atmosphere versus 
the white dwarf's effective temperature is shown in Figure [T2j We see there is no difference 
in the mass in the convective zone, M(Z), between stars with and without an infrared excess. 
An uncertainty is that the mass plotted in Figure [11] are only lower bounds since we do not 
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know how long the accretion has been going on. Also, there is no reason to think that the 
as-yet-to be accreted mass in the disk is a function of the star's temperature. 



5.3. Refractory-rich Parent Bodies? 



Jura fc Xul (120121 ) argue that there are at least two populations of parent bodies that are 



accreted onto a white dwarf: water-rich comet-like objects and dry rocky objects. Here we 
suggest that the compositional variations of the accreted materials might be important in the 
formation of a dust disk. A comparison between calcium, magnesium and iron accretion rates 
is plotted in Figure IT3l log n(Ca)/n(Mg) and log n(Fe)/n(Mg) versus the star's temperature 
are shown in order to correct for different settling rates. These values correspond to the 
intrinsic abundance ratio in the parent body assuming a steady state. 

We see in the upper panel of Figure d3] that most sta rs with an infrared excess have 
relatively high abundance of Ca, a highly refractory element ( lAUegre et al.ll200ll ); the average 
[Ca]/[Mg] in these stars is -0.71 while it is -1.21 in bulk Earth. In cont rast, when measure d 
in a different sample of DZs than studied here, [Ca]/[Mg] equals -1.42 (IKoester et al.ll201ll ). 
In the lower panel of Figure [131 we see Fe abundance almost stays the sam e. Since iron is on e 
of the dominant elements in bulk Earth as well as white dwarf pollution ( iKlein et al.ll2010h . 
there is no reason to think its abundance should vary a lot. There is additional evidence 
supporting that cool white dwarfs are accreting more volatile-rich materials. Sodium, which 
is an important volati le, has a nearly solar abundance in 28 ex treme DZs, much higher than 
the bulk Earth value (IKoester et alihoill : lAUegre et allboOll ). 



6. ORIGIN OF HEAVY ELEMENTS IN STARS WITHOUT AN 

INFRARED EXCESS 



We see most polluted white dwarfs do not possess a dust disk. So where do the heavy 
elements come from? Interstellar accretion was traditionally considered to be the origin of 



this material (jDupuis et al.l 119931 ) but this model has faced many challenges ( iFarihi et al. 



2010al ; IKoester et al.l 120111 ). We discuss two avenues of accretion that can pollute a white 
dwarf's atmosphere without producing a dust disk. 

Accretion from a distant reservoir. It has been suggested that some comets go through 
natural fragmentation due to fast rot ation (Drahus et al.ll201ll) and in this way contribute to 
the zodiacal light in our solar system ( Nesvorny et al.l 2010 ). A similar process might happen 
around a white dwarf given that there might be a large reservoir of comets. Eventually some 
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of the dust particles will d rift inward due to P-R drag and the amount of infrared flux 
produced at frequency v is (I Jura et al.l 120071 ). 



1 



-In 



Ri 
R~ f 



M(Z)c 2 



/47rL> 2 



(3) 



where Ri, R/ is the initial and final distance of the accretion and Rj/R/ is taken to be 1000, 
M(Z) is the total accretion rate, c is the speed of light and D is the distance between the 
Sun and the white dwarf. We compare the flux predicted from Equation (J3]) with our data 
for Cycle 7 targets in Table [3 For the heavily polluted Cycle 7 targets, we see that the 
calculated flux is always much higher than the o bserved value. How ever, there are other cool 
white dwarfs with an accretion rate ~ 10 6 g s _1 (IFarihi et al.ll2009l ) where this process might 
be important. 

Accretion from an orbiting gaseous disk. A gaseous disk that does not emit too much 
in th e infra red can be formed if the grains mutually annihilate as described in the model of 
Jural (120081 ) or the objects sublimate as they get close to the white dwarf. This is likely to 
be tr ue for DAZs warmer than 11,000 K because their settling time is only days (IKoester 
20091 ). Without the continuous feeding of material, the elements would have quickly settled. 



7. CONCLUSIONS 

We find two new stars with an infrared excess, J2209+1223 and WD 0843+516, appar- 
ently the hottest white dwarf with a close-in disk. We fail to find any warm dust around 
heavily polluted DZs, all of which have T* < 10,000 K. We raise the possibility that there 
might be some correlation between the occurrence of dust disk and the volatility of the ac- 
creted material. The best model to explain the source of pollution in white dwarfs without 
an infrared excess has yet to be established. 

We thank J. Farihi for useful email exchanges on Spitzer data, reduction and D. Koester 
for the computation of white dwarf model atmospheres. S. Xu thanks H. Meng and D. 
Rodriguez for helpful suggestions and discussions. This work has been supported by the 
National Science Foundation. 



A. APPENDIX 



In this section, we extend the model of iRafikovl (j2011al ) for P-R drag on the entire disk 
to account for the observed mass accretion rates. Our modification of his model leads to a 
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better agreement with the data. 



Assuming a steady state, the accretion rate provided by P-R drag is (lRafikovll2011af ). 



^^{j" inner) 



160 r rl oTj 

3 Tin,np,T C 



(Al) 



where a is the Stefan Boltzmann constant, c is the speed of light, and <p r is an efficiency 
coefficient. 



Using Equation (jSJ), iRafikovl feOllal ) derived an accretion rate 



32cj) r ( r*T*T s 
— - — o 



1.80 x 10 s 



T* V _! 
J 9 s 



15000fsT 



(A2) 



In the last step we take r*=0.013 R and <f> r = 1. 



As discussed in section |4"T2"| our inner boundary is described by Equation ([T]); so we get 
a different value of accretion rate, 



M{r ir , 



3tt 



1/3 



I6<p r a {r*\ 2 



8.14 x 10 s 



3 v c 
15000fsT 



(TsT 2 } 4/3 



\ 8/3 

^ g*- 1 



(A3) 



A comparison between these two accretion rates are shown in Figure UM and we see they 
differ as much ClS db factor of 5. 

Observationally, we can also derive the accretion rate from the atmospheric pollution. 
Since we are trying to make a comparison between models that vary by a factor of 5, we 
need accurate abundances. Therefore, we only consider stars t hat have measurem ents for at 
least three of the four major elements in white dwarf pollution (IKlein et al.ll2010l ): O, Mg, Si 
and Fe. In Figure dH we see the accretion rates derived from our model and spectroscopical 
analysis agree well for DAZs. Since the settling time for heavy elements is relatively short 
in the se stars, the p resence of a dusty disk indicates the accretion should be in a steady 
state (iKoesterl 120091 ). which satisfies the assumption of the calculation above. For DBZs, 
the settling time is comparatively longer and we cannot determine if the system is in the 
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build-up phase or steady state. The accretion rate driven by PR drag gives a good lower 
bound in these stars and there might be other mechanisms on top of that, as discussed in 



Rafikovl fl2011bf ). 



REFERENCES 

Allegre, C, Manhes, G., & Lewin, E. 2001, Earth and Planetary Science Letters, 185, 49 

Bohlin, R. C, Gordon, K. D., Rieke, G. H., Ardila, D., Carey, S., Deustua, S., Engelbracht, 
C, Ferguson, H. C, Flanagan, K., Kalirai, J., Meixner, M., Noriega-Crespo, A., Su, 
K. Y. L., & Tremblay, P. 2011, ArXiv e-prints 

Bonsor, A., Mustill, A. J., & Wyatt, M. C. 2011, MNRAS, 594 

Bonsor, A. & Wyatt, M. 2010, MNRAS, 409, 1631 

Brinkworth, C. S., Gansicke, B. T., Marsh, T. R., Hoard, D. W., & Tappert, C. 2009, ApJ, 
696, 1402 

Chu, Y., Gruendl, R. A., Bilikova, J., Riddle, A., & Su, K. 2010, ArXiv e-prints 

Chu, Y.-H., Su, K. Y. L., Bilikova, J., Gruendl, R. A., De Marco, O., Guerrero, M. A., 
Updike, A. C, Volk, K., & Rauch, T. 2011, ArXiv e-prints 

Debes, J. H., Hoard, D. W., Kilic, M., Wachter, S., Leisawitz, D. T., Cohen, M., Kirkpatrick, 
J. D., & Griffith, R. L. 2011, ApJ, 729, 4 

Debes, J. H. & Sigurdsson, S. 2002, ApJ, 572, 556 

Dong, R., Wang, Y., Lin, D. N. C, & Liu, X.-W. 2010, ApJ, 715, 1036 

Drahus, M., Jewitt, D., Guilbert-Lepoutre, A., Waniak, W., Hoge, J., Lis, D. C, Yoshida, 
H., Peng, R., & Sievers, A. 2011, ApJ, 734, L4+ 

Dufour, P., Bergeron, P., Liebert, J., Harris, H. C, Knapp, G. R., Anderson, S. F., Hall, 
P. B., Strauss, M. A., Collinge, M. J., & Edwards, M. C. 2007, ApJ, 663, 1291 

Dufour, P., Bergeron, P., Schmidt, G. D., Liebert, J., Harris, H. C, Knapp, G. R., Anderson, 
S. F., & Schneider, D. P. 2006, ApJ, 651, 1112 

Dufour, P., Kilic, M., Fontaine, G., Bergeron, P., Lachapelle, F., Kleinman, S. J., & Leggett, 
S. K. 2010, ApJ, 719, 803 



- 14 - 



Dupuis, J., Fontaine, G., Pelletier, C, & Wesemael, F. 1993, ApJS, 84, 73 

Eisenstein, D. J., Liebert, J., Harris, H. C, Kleinman, S. J., Nitta, A., Silvestri, N., An- 
derson, S. A., Barentine, J. C, Brewington, H. J., Brinkmann, J., Harvanek, M., 
Krzesinski, J., Neilsen, Jr., E. H., Long, D., Schneider, D. P., & Snedden, S. A. 2006, 
ApJS, 167, 40 

Farihi, J., Barstow, M. A., Redfield, S., Dufour, P., & Hambly, N. C. 2010a, MNRAS, 404, 
2123 

Farihi, J., Becklin, E. E., & Zuckerman, B. 2008a, ApJ, 681, 1470 

Farihi, J., Brinkworth, C. S., Gansicke, B. T., Marsh, T. R., Girven, J., Hoard, D. W., Klein, 
B., & Koester, D. 2011, ApJ, 728, L8 

Farihi, J., Jura, M., Lee, J., & Zuckerman, B. 2010b, ApJ, 714, 1386 

Farihi, J., Jura, M., & Zuckerman, B. 2009, ApJ, 694, 805 

Farihi, J., Zuckerman, B., & Becklin, E. E. 2008b, ApJ, 674, 431 

Fazio et al. 2004, ApJS, 154, 10 

Girven, J., Gaensicke, B. T., Steeghs, D., & Koester, D. 2011, ArXiv e-prints 

Hoard, D. W., Wachter, S., Sturch, L. K., Widhalm, A. M., Weiler, K. P., Pretorius, M. L., 
Wellhouse, J. W., & Gibiansky, M. 2007, A J, 134, 26 

Holberg, J. B. & Bergeron, P. 2006, AJ, 132, 1221 

Jura, M. 2003, ApJ, 584, L91 

-. 2008, AJ, 135, 1785 

Jura, M., Farihi, J., & Zuckerman, B. 2007, ApJ, 663, 1285 
Jura, M. & Xu, S. 2012, AJ, submitted 

Kilic, M., Farihi, J., Nitta, A., & Leggett, S. K. 2008, AJ, 136, 111 

Kilic, M., Kowalski, P. M., Reach, W. T., & von Hippel, T. 2009, ApJ, 696, 2094 

Kilic, M., Munn, J. A., Harris, H. C, Liebert, J., von Hippel, T., Williams, K. A., Metcalfe, 
T. S., Winget, D. E., & Levine, S. E. 2006, AJ, 131, 582 



- 15 - 



Klein, B., Jura, M., Koester, D., & Zuckerman, B. 2011, ApJ, in press 

Klein, B., Jura, M., Koester, D., Zuckerman, B., & Melis, C. 2010, ApJ, 709, 950 

Koester, D. 2009, A&A, 498, 517 

Koester, D., Girven, J., Gansicke, B. T., & Dufour, P. 2011, A&A, 530, A114+ 

Koester, D., Voss, B., Napiwotzki, R., Christlieb, N., Homeier, D., Lisker, T., Reimers, D., 
& Heber, U. 2009, A&A, 505, 441 

Liebert, J., Bergeron, P., & Holberg, J. B. 2005, ApJS, 156, 47 

Martin, D. C, Fanson, J., Schiminovich, D., Morrissey, P., Friedman, P. G., Barlow, T. A., 
Conrow, T., Grange, R., Jelinsky, P. N., Milliard, B., Siegmund, O. H. W., Bianchi, 
L., Byun, Y., Donas, J., Forster, K., Heckman, T. M., Lee, Y., Madore, B. F., Malina, 
R. F., Neff, S. G., Rich, R. M., Small, T., Surber, F., Szalay, A. S., Welsh, B., & 
Wyder, T. K. 2005, ApJ, 619, LI 

Melis, C, Farihi, J., Dufour, P., Zuckerman, B., Burgasser, A. J., Bergeron, P., Bochanski, 
J., & Simcoe, R. 2011, ApJ, 732, 90 

Melis, C, Jura, M., Albert, L., Klein, B., & Zuckerman, B. 2010, ApJ, 722, 1078 

Nesvorny, D., Jenniskens, P., Levison, H. F., Bottke, W. F., Vokrouhlicky, D., & Gounelle, 
M. 2010, ApJ, 713, 816 

Rafikov, R. R. 2011a, ApJ, 732, L3+ 

-. 2011b, MNRAS, L287+ 

Reach, W. T., Megeath, S. T., Cohen, M., Hora, J., Carey, S., Surace, J., Willner, S. P., 
Barmby, P., Wilson, G., Glaccum, W., Lowrance, P., Marengo, M., & Fazio, G. G. 
2005, PASP, 117, 978 

Reid, I. N., Liebert, J., & Schmidt, G. D. 2001, ApJ, 550, L61 

Steele, P. R., Burleigh, M. R., Dobbie, P. D., Jameson, R. F., Barstow, M. A., & Satterth- 
waite, R. P. 2011, ArXiv e-prints 

Vennes, S., Kawka, A., & Nemeth, P. 2011, MNRAS, 413, 2545 

von Hippel, T., Kuchner, M. J., Kilic, M., Mullally, F., & Reach, W. T. 2007, ApJ, 662, 544 
Werner et al. 2004, ApJS, 154, 1 



- 16 - 

Zacharias, N., Monet, D. G., Levine, S. E., Urban, S. E., Gaume, R., & Wycoff, G. L. 2005, 
VizieR Online Data Catalog, 1297, 

Zuckerman, B. & Becklin, E. E. 1987, Nature, 330, 138 

Zuckerman, B., Koester, D., Reid, I. N., & Hiinsch, M. 2003, ApJ, 596, 477 

Zuckerman, B., Melis, G, Klein, B., Koester, D., & Jura, M. 2010, ApJ, 722, 725 



This preprint was prepared with the A AS IATgX macros v5.2. 



-17- 



Table 1. Warm Dust Disks and Candidates Discovered Since 2010 



WD 


Name 


SpT 


T* 


V 


Discovery 


Discovery 


Ref 








(K) 


(mag) 


Year 


Telescope 




0106-328 


HE 0106-3253 


DAZ 


15,700 


15.50 


2010 


Spitzer 


1 


0307+077 


HS 0307+0746 


DAZ 


10,200 


16.40 


2010 


Spitzer 


1 


0435+410 


GD 61 


DBZ 


17,280 


14.80 


2011 


Spitzer 


2 




J0738+1835 


DBZ 


13,600 


17.82 d 


2010 


Germini/NIRI 


3 


0843+516 


PG 0843+517 


DA 


23,900 


16.15 


2011 


Spitzer 


4 


1225-079 a 


PG 1225-079 


DZAB 


10,500 


14.80 


2010 


Spitzer 


1 




J2209+1223 


DBZ 


17,300 


17.43 d 


2011 


Spitzer 


4 




GALEX J1931+0117 


DAZ 


20,890 b 


14.20 


2011 


WISE 


5 


2221-165 


HE 2221-1630 


DAZ 


10,100 


16.10 


2010 


Spitzer 


1 



SDSS 0753+2447 c 


DA 


13,400 


19.21 d 


2011 


UKIRT 


6,7 


SDSS 0959-0200 c 


DA 


12,000 


18.34 d 


2011 


UKIRT 


6 


SDSS 1221+1245 c 


DA 


12,000 


18.39 d 


2011 


UKIRT 


6 


J1320+0018 c 


DA 


19,600 


17.30 


2011 


UKIRT 


7 


J1557+0916 c 


DA 


22,000 


18.17 d 


2011 


UKIRT 


7 



a WD 1225-079 is reported in lFarihi et all (j2010bl ) to have a 4cr excess in the IRAC 7.9 /mi band. 



Melis et all ( 1201 ll ) derived a stellar temperature of 23,470 K for this star. 



c These are debris disk candidates that needs to be confirmed with more infrared data as the 
UKIRT observation only extends to the K band. 

d This is SDSS r magnitude. 



Note. - This table is supplementary to Table 1 in iFarihi et al.l ( 120091 ). which lists 14 white 
dwarfs with an infrared excess discovered at that time. 



Note. - 



GD 303, which might have a marginal excess at la level at the IRAC 4 band ( IFarihi et al 
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2010bl ). is excluded from this table. 
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Table 2. Cycle 7 White Dwarf Targets 







T* 


SDSS z 






WD 


Name 


(K) 


(mag) 


T /^i 1 / Ttt In 

[Ca]/[He] a 


Ref 




T003fif)1 37-1 1 1 913 8 


7 980 


17 1 8 

1 1 .10 


-Q 9fi 


1 

1 


U^l U Uc/J 


T09183fi RQ-0Q1Q44 8 


Q FifiD 


1 7 84 

J. 1 .01 


-1 n fis 


i 

1 




T0Q3Q49 9Q+555048 7 


o,uou 


1 1 . ±o 


-8 51 


1 

1 






R ^70 


1 7 7^ 


1 n 97 


1 

1 


1035-003 


J103809.19-003622.5 


8,200 c 


17.30 


-7.70 


2 


1212-023 


J121456.39-023402.8 


6,000 


17.50 




3 


1244+498 


J124703.28+493423.6 


16,800 


17.21 




4 




J125752. 77+425255.1 


16,800 


17.94 




4 




J130905.26+491359.7 


8,620 


17.67 


-10.16 


1 




G 199-63 


8,900 


17.15 


-9.39 


1 




LP 219-80 


6,770 


16.91 


-11.35 


1 




J220934.84+122336.5 


17,300 


17.91 


-6.16 d 


4 




J222802.05+120733.3 


6,760 


16.36 


-9.96 


1 


2222+683 


G 241-6 


15,300 


16.07 


-7.25 


5 



a Abundances are expressed as [X]/[Y] = log [n(X)/n(Y)]; n(X) is the num- 
ber abundance. 



b This star is reported to be 11,500 K in iKilic et all (120081 ) 



This star is reported to be 6,770 K and [Ca]/[He]= -9.44 in iDufour et al. 



(120071 ) 



d This value is derived from measuring the equivalent width of Ca II 3933. 7A 
from the SDSS spectra and extrapolate from GD 61, a heavily polluted white 
dwarf with similar effective temperature and surface gravity. 
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Table 3. Archived White Dwarf Targets 





rp 

J-eff 


SDSS z 




WD 


(K) 


(mag) 


Ref 


0816+297 


16 700 

-L W ■ 1 WW 


1 6 53 

X W > WW 


1 


Ww ± U \^ O W (J 


1 8 700 

-LtJ j 1 WW 


1 6 38 


1 


0843+516 


23,900 


16.86 


1 


0937+505 


35,900 


16.86 


1 


1017+125 


21,400 


16.55 


2 


1109+244 


37,800 


16.64 


3 


1120+439 


27,200 


16.28 


3 


1133+293 


23,000 


15.69 


3 


1214+267 


65,700 


16.54 


3 


1216+036 


14,400 


16.66 


2 


1257+032 


17,600 


16.39 


2 


1507+021 


20,200 


17.23 


2 


1553+353 


25,600 


15.59 


3 


1559+128 


29,200 


17.74 


3 


1620+513 


20,900 


16.61 


1 


2120+054 


36,200 


16.93 


1 
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Table 4. IRAC Fluxes for Cycle 7 Targets 





F3.6/*m 


F4.5/«n 


Name 


/xJy 


/xJy 


WD 0033-114 


61+4 


41+3 


WD 0216-095 


28+2 


17+2 


WD 0936+560 


60+4 


37+3 


J0951+4033 


34+3 


22+2 


WD 1035-003 


53+3 


36+3 


WD 1212-023 


85+5 


60+4 


WD 1244+498 


38+3 


23+2 


J1257+4252 


19+2 


13+2 


J1309+4913 


38+3 


24+2 


G 199-63 


57+4 


37+3 


LP 219-80 


104+7 


71+5 


J2209+1223 


88+5 


86+5 


J2228+1207 


156+9 


104+6 


WD 2222+683 


95+5 


64+4 
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Table 5. IRAC Fluxes for Archived DAs 







4.5/zm 




r 7.9fim 


WD 


u,Tv 


//.Tv 


//.Tv 


//.Tv 


0816+297 


76+5 


51+4 


28 a 


25 a 


0819+363 


97+5 


50+3 


66+8 


55+10 


0843+516 


140+12 


134+9 


98+10 


154+13 


0937+505 


43+3 


27+3 


22+7 


27 a 


1017+125 


71+6 


47+6 


35 a 


24 a 


1109+244 


57+5 


36+4 


37 a 


29 a 


1120+439 


83+5 


47+4 


65+11 


54 a 


1133+293 


147+9 


90+6 


61+13 


53+16 


1214+267 


54+4 


35+3 


32 a 


19 a 


1216+036 


74+4 


45+3 


34 a 


36 a 


1257+032 


85+5 


51+4 


34+11 


33 a 


1507+021 


40+3 


25+3 


32 a 


26 a 


1553+353 


167+9 


102+6 


60+7 


38+9 


1559+128 


21+2 


12+2 


38 a 


21 a 


1620+513 


69+4 


43+3 


23 a 


26 a 


2120+054 


45+3 


27+3 


39 a 


23 a 



3cr upper limit. 



Note. The WISE photometry for WD 

0819+363 and WD 1559+128 might be somewhat 
misleading as the IRAC images show a background 
source within 5'.'0, which cannot be resolved by 
WISE and its measurement indicates a substantial 
infrared excess for these two stars. Confusion is 
also likely t o be a p roblem for GALEX 1931+0117, 



Melis et al.l ( 1201 ll ) found a significantly lower L- 
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band flux with g round based photometry than 



Debes et all (120 111 ) did by using WISE. 



Table 6. Disk Model Parameters 



Name 




r*/D a 


T 

1 inner 


Touter 


dinner 


^outer 


cos i 




(K) 


(io- 12 ) 


(K) 


(K) 


(R*) 


(R*) 




WD 0843+516 


23,900 


3.9 


1,260 


640 


30 


75 


0.17 


J2209+1223 


15,000 


1.6 


1,340 


590 


15 


45 


0.08 






1.6 


910 


710 


25 


35 


0.29 



This is the ratio of the stellar radius and its distance to the Sun. 



Table 7. Infrared Flux Produced by Accretion from Cometary Dust 



Name 


D 


log M(Z) 


F/U/wn^Jy) 






(PC) 


(g/s) 


predicted 11 


observed 


WD 0033-114 


55 


8.66 


6,300 


41+3 


WD 0216-095 


97 


7.34 


120 


17+2 


WD 0936+560 


65 


9.41 


25,000 


38+2 


J0951+4033 


84 


7.69 


310 


22+2 


WD 1035-003 


49 


9.65 


77,000 


37+2 


J1309+4913 


83 


7.76 


370 


24+2 


G 199-63 


68 


8.53 


3,100 


37+2 


LP 219-80 


45 


6.61 


160 


71+4 


J2228+1207 


34 


7.96 


3,200 


105+6 



a These values are the predicted total flux, which is the sum of 
the photospheric flux and the excess due to accretion from the 
cometary dust. 



Note. - Thi s table contains all the Cycle 7 stars with re- 
ported distance (IDufour et al.ll2006l ). 
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0.059 0.063 0.06B 0.072 0.077 0.0B2 0.086 0.091 0.095 



Fig. 1. — IRAC mosaics of J2209+1223 in 3.6 /j,m (the upper panel) and 4.5 /zm band (the 
lower panel) in the plate scale of 1'.'2 pixel -1 ; the left columns are the original data and the 
right columns are the residue after PRF fitting and the white dwarf is subtracted. North is 
up and East is left; the field of view is 36'.'0 by 36'.'0. The red circle is centered on J2209+1223 
with a radius of 3 pixels and the nearby star is 6'.'2 away at position angle 333°. 
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0.038 0.Q4B 0.054 0.062 0.O7 0.078 0.08B 0.094 0.1 




0.0351 0.0412 0.0473 0.0534 0.0595 0.0856 0.0717 0.0778 0.0839 



Fig. 2. — The same as Figure [T] except for G 241-6. We see this is a very complicated field 
and G 241-6 is heavily blended with an object with an angular separation of only 2'.'3 at 
position angle 170° and another galaxy, J222334. 23+683723. 8 with a separation of 5'.'5 and 
position angle 110°. 
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Fig. 3.— SED for Cycle 7 targets, including data from SDSS, 2MASS and IRAC with 2a 
error bars. The dashed line is a simple blackbody fit to the photospheric flux. 
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Fig. 4. — The same as Figure [3j 
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Fig. 5. — The same as Figure [3] except for G 241-6. 
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Fig. 6.— SED for archived targets, including data from SDSS, 2MASS, IRAC and WISE 
data when available with 1o error bars. The asterisks denote upper limit. The dashed line 
is a simple blackbody fit to the photospheric flux. 
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Fig. 7. — The same as Figure EJ 



-32 - 



WD 1559+128 
29200 K 



z i r g u 



-1 -0.8 -0.6 -0.4 -0.2 

log (nrrf 1 ) 



0.2 0.4 0.6 



WD 1620+513 
20900 K 



z i r g u 



K H J 



W2 W1 



jt ■ » - — 



-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 



log X 1 (nm 1 ) 



3500 
3000 
2500 
2000 
1500 



WD 2120+054 


I I 1 

z i r g u 


36200 K 




K H J 


/ J 

• t 






*. , * .^-w-.* ----- 





-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 O.f 

log X -1 (urn -1 ) 



Fig. 8. — The same as Figure EJ 



-33 - 




r 1 (urn -1 ) 




r 1 (urn -1 ) 

Fig. 9. — SED for J2209+1223. The green line displays the photospheric flux from the star. 
The black lines and dashed pink lines represent the two disk models listed in Table |6j the 
flux from the disk and its total flux. The upper panel shows the fit to the entire model while 
the lower panel shows the fit to the disk; 2cr error bars are displayed. 
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Fig. 10. — SED for WD 0843+516. The green line displays the model atmosphere (D. 
Koester, private communication), the blue line the model disk and the black line the sum of 
both. The upper panel shows the entire fit to the model atmosphere while the lower figure 
only shows the disk portion. The pink dots are taken from GALEX, SDSS, 2MASS, green 
dots from WISE and red dots from IRAC, respectively; 2a error bars are displayed. 
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Fig. 11. — Accretion rates at different stellar temperatures. The data include stars with well 



deter mined m ajor element abundance , newly identified heavily polluted DZs (IKoester et al. 



201ll ). DAZs (jZuckerman et al.ll2003l ) and Cycle 7 stars reported in this paper . The dash- 
dotted line represents an accretion rate of 3 x 10 8 g s _1 , a threshold proposed by lFarihi et al. 



( 120091 ): over 50% of single white dwarfs with accretion rate at least this high display an 
infrared excess from ~1000 K dust. 
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Fig. 12. — Mass of the pollution in the convective zone versus temperature for all the DBZs 
and DZs in Figure [TTJ The black dashed line denotes the mass in the mixing zo ne of GD 61, 



1.6 x 10 g, which is the minimum mass found so far to have an infrared excess (IFarihi et al. 



201ll ). We see most of the stars are accreting from parent bodies more massive than this 
value. 
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Fig. 13. — Compositional difference for all the stars in Figure [TT] that have Ca, Mg and Fe 
detections. The relative abundance has been corrected for settling assuming a steady state. 
The dashed line is the value for bulk Earth. We see the stars with infrared excess usually 
have a high n(Ca)/n(Mg). 
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Fig. 14. — Mass accretion rate onto white dwarfs with an infrared excess. The s olid line 



shows t he value predicted by ou r model while the dashed line is th e value derived inlRafikov 



f|2011al ). The DA Zs are G29-38 (IKoesterll2009h and GALEX 19 31 flVennes et a 



DBZ s are GD 40 fclein et aliboioh. G D 61 iFarihi et aDboilh . SPSS J0738 JDufour et aL 



2011). And 



2010|), PGJj2^ TOjK lem et l201lh . Ton 345 (D. Koester, private comminucation) and 



GD 362 (Koester 20091 ). All the stars in this figure have well determined major element 
abundance. 



